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£j , ABSTRACT 

We report on observations of the X-ray pulsar IGR J16320-4751 (a.k.a. AX J163 1.9-4752) 
performed simultaneously with INTEGRAL and XMM-Newton. We refine the source posi- 
tion and identify the most likely infrared counterpart. Our simultaneous coverage allows us 
to confirm the presence of X-ray pulsations at ~ 1300 s, that we detect above 20 keV with 
r> | INTEGRAL for the first time. The pulse fraction is consistent with being constant with en- 

^ ■ er gy> which is compatible with a model of polar accretion by a pulsar. We study the spectral 

properties of IGR J16320-4751 during two major periods occurring during the simultaneous 
coverage with both satellites, namely a flare and a non-flare period. We detect the presence 
of a narrow 6.4 keV iron fine in both periods. The presence of such a feature is typical of 
supergiant wind accretors such as Vela X-l or GX 301-2. We inspect the spectral variations 
with respect to the pulse phase during the non-flare period, and show that the pulse is solely 
due to variations of the X-ray flux emitted by the source and not to variations of the spec- 
tral parameters. Our results are therefore compatible with the source being a pulsar in a High 
Mass X-ray Binary. We detect a soft excess appearing in the spectra as a blackbody with a 
temperature of ~0.07 keV. We discuss the origin of the X-ray emission in IGR J 16320-4751: 
while the hard X-rays are likely the result of Compton emission produced in the close vicinity 
of the pulsar, based on energy argument we suggest that the soft excess is likely the emission 
by a collisionally energised cloud in which the compact object is embedded. 

Key words: accretion - x-rays: binaries - stars: neutron - stars: pulsar:general - stars: individ- 
ual IGR J 1 6320-475 1 , AX J 1 63 1 .9-4752 
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1 INTRODUCTION 

The INTErnational Gamma-Ray Astrophysics Laboratory (INTE- 
GRAL) was launched on October 17, 2002 (Winkler et al. 2003). 
Since then, through regular scans of our Galaxy and guest ob- 



servers' observations, about 75 new sources 1 have been discov- 
ered mainly with the IBIS telescope (Ubertini et al. 2003), thanks 

1 An updated list of all INTEGRAL sources can be found at 
http://isdc.unige.ch/~rodrigue/html/igrsources.html 
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to its low energy camera ISGRI (Lebrun et al. 2003). Because of 
its energy range (from 15 keV to ~ 1 MeV), its high angular 
resolution (12 arcmin), good positional accuracy (down to ~ 0.5 
arcmin for bright sources), and its unprecedented sensitivity be- 
tween 20 and 200 keV, IBIS/ISGRI has helped to answer the ques- 
tion of the origin of the hard X-ray background in the Galaxy 
iLebrun et al. 20041) . These capacities have also allowed to discover 
many peculiar X-ray binaries characterized by a huge equivalent 
absorption column de nsity (Nyt), as high as a few times 10 24 cm - 2 
in IGR J 1 63 1 8-4848 iMatt & Guainazzi 2003tlWalter et al. 2o"o3h . 
the first source discovered by INTEGRAL. Due to the high absorp- 
tion, most of these sources were not detected during previous soft 
X-ray scans of the Galaxy (see e.g. Kuulkers 2005 for a review). 

IGR J 16320-4751 w as detected on Feb. 1, 2003, with IS- 
GRI iTomsick et al. 20031) as a hard X-ray source. The source was 
observed to vary significantly in the 15-40 keV energy range on 
time scales of ~ 1000 s, and was sometimes detected above 60 
keV iTomsick et al. 20031 iFoschini et al. 20041) . Inspection of the 
X-ray archives revealed that IGR 116320-4751 is the hard X-ray 
counterpa rt to AX J1631.9-47 52, observed with ASCA in 1994 
and 1997 JSugizaki et al. 200 ll) . Analysis of archival BeppoSAX- 
WFC data showed that this source was persistent for at least 
8 years (in't Zand et al. 2003). Soon after the discovery of 
IGR J16320-4751 by INTEGRAL, an XMM-Newton Target of Op- 
portunity was triggered. Thi s allowed us to obtain the most accu- 
rate X-ray position to date jRodriguez et al. 20031) . which in par- 
ticula r led to the identification of two possible infrared counter- 
parts iTomsick et al. 2003llRodriguez et al. 20031) (hereafter source 
1 and 2). From this analysis, we suggested that IGR J16320-4751 
is proba bly a High Mass X-ra y Binary (HMXB) hosting a neu- 
tron star iRodriguez et al. 20031) . This last point has been reinforced 
since the discovery of X-ray pulsatio ns from this source in both 
XMM-Newton and ASCA observations iLutovinov et al. 2005T) with 
a pulse period of about 1300 s. Aharonian et al. (2005) reporting 
results of the survey of the inner regions of the Galaxy at very high 
energy with HESS, found a new source, HESS J1632-478, at a po- 
sition coincident with IGR J 16320-4751, but the authors suggest 
that this could be simply a chance coincidence. Very recently, Cor- 
bet et al. (2005) reported the discovery of strong modulations of 
the X-ray flux seen with Swift, that they interpreted as revealing 
the orbital period of the system. The period of 8.96 days that they 
find is compatible with the system containing an early type super- 
giant (Corbet et al. 2005). 

We report here observations of IGR J 16320-4751 performed si- 
multaneously with INTEGRAL and XMM-Newton in August 2004. 
The sequence of observations is introduced in the following section 
together with technical information concerning the data reduction. 
Section 3 of this paper describes the results which are commented 
and discussed in the last part of the paper. 



2 OBSERVATIONS AND DATA REDUCTION 

The journal of the observations can be found in TableQ The INTE- 
GRAL observation results from an amalgamation of the observation 
of IGR J 163 18-4848 (PI Kuulkers, programme # 0220007) with 
that of IGR J16320-4751 (PI Foschini, programme # 0220014). 
IGR J 16320-4751, however, remained the on-axis target for this 
observation. 



2.1 INTEGRAL Data Reduction 

Our INTEGRAL observatio n was performed with the so-called 
hexagonal dithering pattern ICourvoisier et al. 20031) . which con- 
sists of a sequence of 7 pointings (called science windows, here- 
after: SCW) following a hexagonal dithering on and around the po- 
sition of the source. Being the on-axis target, IGR J16320-4751 is 
always in the fully coded field of view of the IBIS and SPI instru- 
ments, where the instrumental response is optimal. The INTEGRAL 
data were reduced using the Off-line Scientific Analysis (OSA) v 
5.0, with specificities for each instrument described below. 

The data from IBIS/ISGRI were first processed until produc- 
tion of images in the 20-40 and 40-80 keV energy ranges, with the 
aim of identifying the most active sources in the field. From the re- 
sults of this step, we produced a catalogue of active sources which 
was given as an input for a second run producing images in the 
20-60 and 60-200 keV energy ranges. In the latter we forced the 
extraction/cleaning of each of the catalogue sources in order to ob- 
tain the most reliable results for IGR J 16320-475 1 (see Goldwurm 
et al. 2003 for a detailed description of the IBIS software). In order 
to check for the presence of the 1300 s X-ray pulsations in the IBIS 
range, we also produced light curves with a time bin of 250 s (with 
the 0SA5 . Q/ii_lc_extract v2 . 4 . 3 module) in the same energy 
ranges. 

We produced spectra using two alternative methods. The first 
method is to use the standard spectral extraction from the OSA 
pipeline. The second method uses the individual images produced 
in 20 energy bins, to estimate the source count rate and build 
the spectrum as explained in Rodriguez et al. (2005). This sec- 
ond method can be used to cross check the results obtained with 
the standard spectral extraction. Comparison of the spectra ob- 
tained with the two methods showed no significant differences, 
we therefore used the spectra obtained with the standard proce- 
dure in our spectral fits. We used the standard response matri- 
ces provided with OSA 5.8, i.e. isgr_arf_rsp_0010 . fits and 
isgr_rmf_grp_0Q16 .fits, the latter rebinned to 20 spectral chan- 
nels. 

The source is not spontaneously detected by the IEM-X de- 
tector. We tentatively forced the extraction of science products at 
the position of IGR J16320-4751, but a rapid comparison with 
the XMM-Newton data showed that the JEM-X products were 
very likely contaminated by the nearby black hole candidate 4U 
1630-47 that was in a bright soft state outburst at the time of the 
observation (Tomsick et al. 2005). These data were not considered 
further. We did not use the data from the spectrometer SPI, because 
the 2.5° angular resolution did not allow to discriminate the emis- 
sion of IGR J16320-475 1 from that of e.g. 4U 1630-47 (the angu- 
lar separation between the 2 sources is ~ 40'), or IGR J163 18-4848 
(~ 55'). 



2.2 XMM-Newton Data Reduction 

The XMM-Newton data were reduced with the Science 
Analysis System (SAS) v6.1.0. Event lists for EP IC MOS 
iTurner et al. 200 ll) and EPIC PN IStriider et al. 200 ll) cameras 
were obtained after processing of the ODF data with emchain 
and epchain. During the processing, the data were screened 
by rejecting periods of high background, and by filtering the 
bad pixels. The EPIC MOS were both operating in timing mode 
allowing one to obtain light curves from the central chip with 1.5 
ms resolution. The EPIC PN was operating in full frame mode, 
allowing light curves with a time resolution as high as 73.4 ms 
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Table 1. Journal of the observations analysed in this paper. 



Satellite 


Revolution 




Date obs. 




Prop. Id 


total duration 




# 


Start day 


Stop day 


(MJD-53000) 




(Xl0 3 s) 


INTEGRAL 


226 


2004-Aug-19 


2004-Aug-20 


236.57-237.94 


0220014& 0220007 


~ 118 


XMM-Newton 


860 


2004-Aug-19 


2004-Aug-20 


236.55-237.14 


0201700301 


~ 51 



to be obtained. For the PN camera we extracted the spectra and 
light curves from a 35" radius circle centered on the source, while 
background products were extracted from a 60" circular region 
free of sources (from the same chip). 

Because of the operating mode of the MOS cameras, no 
background estimate can be obtained from the central chip where 
the source lies. A "quick look" to the lateral MOS chips has 
shown that the background remained negligible during the whole 
observation. The latter could therefore be neglected in the analysis 
of the MOS data. MOS light curves were hence extracted from 
the central chip of both cameras. For the three EPIC cameras 
light curves were extracted in different energy ranges (2-12 keV, 
0.6-2 keV, 2-5 keV, 5-8 keV, 8-12 keV) with the highest time 
resolution achievable. Barycentric correction was applied to all 
light curves. The MOS light curves were then summed using the 
FT00LS lcmath. 

A redistribution response matrix and an ancillary response 
file for the PN spectral products were generated with rmfgen 
and arfgen. T he resulting spectrum was fitted in XSPEC vl 1.3.1 
jAmaud 199(3) simultaneously with the spectrum obtained with 
INTEGRAL/IS GRI . 



Figure 1. 2MASS K s 1.3'x0.9' image of the field around 
IGR J16320-4751. The 2 circles represent respectively the XMM- 
Newton position given in Rodriguez et al. (2003, large circle), and the 
refined position discussed in the text. Infrared sources 1 and 2 are the 
candidate counterparts. 




3 RESULTS 

3.1 Refining of the X-ray position 

Given the long exposure time (Table Q and the high flux from 
the source we tried to obtain a better estimate on the position 
in order to possibly discrimina te between the 2 candidate coun- 
terparts iRodriguez et al. 20031) . For this purpose we used the 
edetect_chain task after having extracted PN images in 5 en- 
ergy ranges. The latter were further rebinned so that an image 
pixel had a physical size of 4.4". The best position obtained with 
this method is RA 72000 =248.0077° and DECy 2 ooo=-47.8742° with 
a nominal uncertainty of ~ 4". This is about 1.9" from the posi- 
tion reported by Rodriguez et al. (2003). To cross check, we re- 
analysed the 2003 XMM-Newton observation with the latest cali- 
bration files available, and using the same edetect.chain in the 
same energy ranges. As explained in Rodriguez et al. (2003), only 
MOS data are available, and the data need to be filtered for soft 
proton flares. The best position we obtained is RAj2ooo=248.0081° 
and DECj2ooo= - 47.8741°. We therefore can refine the X-ray po- 
sition to RAj2ooo = 16 A 32'" 01*9 DEC y2 ooo= -47° 52' 27" with an 
uncertainty of ~ 3" at 90% confidence. This new position suggests 
more strongly that the source labeled 1 in Fig. is the genuine 
counterpart of IGR J 16320-4751, since source 2 is now outside of 
the 90% error box on the X-ray position. 

We compared the X-ray source positions for high significance 
sources (likelihood ^ 200) with infrared sources from the 2MASS 
catalog. We found a match within 0.4" between the XMM-Newton 
source at RA J20 oo=16' ! 32"' 12 ( .306, DEC /2000 =-47 ° 44' 59".32 
and a bright 2MASS source with a K magnitude of 9.3 lo- 
cated at RA y20 oo=16'' 32'" 12*339, DEC 20 oo=-47° 44' 59".45. We 



found a second close match between the XMM-Newton source at 
RA J20 oo=16'' 31"' 351642, DEC, 20 o=-47° 51' 27". 42 and a bright 
2MASS source with a K magnitude of 9.3 located at RAj 20 oo=16'' 
31"' 351634, DEC J20 oo=-47° 51' 27".89. Allowing a 1" error circle 
for the XMM-Newton source and taking into account the surface 
density of XMM-Newton sources and 2MASS sources at appropri- 
ate magnitude limits, we find that the chance probability of occur- 
rence of these X-ray/IR source coincidences is 0.3%. This increases 
our confidence that the astrometry of the X-ray is correct and that 
we have identified the correct IR counterpart to IGR J16320-4751. 



3.2 Timing Analysis 

The light curves of IGR J 16320-4751 obtained in different energy 
ranges with the different instruments are presented in Fig. [5] 
As already reported, IGR J16320-4751 is a variable source on 
various time scales, from seconds-minutes (Rodriguez et al. 2003, 
Fig. to days-months (e.g. Foschini et al. 2004). During our 
observations, the source shows a prominent flare around MJD 
53236.6, visible in both light curves (Fig. |2j- A second similar 
flare is visible in the ISGRI 20-60 keV light curve around MJD 
53237.3. Unfortunately, our XMM-Newton observation does not 
cover this period. 

We used the EPIC cameras to study the pulsations of 
IGR J 16320-4751 already reported in Lutovinov et al. (2005). We 
performed a period search on the MOS and PN light curves using 
the XROHO S tool efsea rch. Because the source is highly absorbed 
iRodriguez et al. 20031) . and in order to improve the detection of 
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Figure 2. Left: The 2-12 keV XMM-Newton/PN light curve of IGR J 1 6320-475 1 . Right : The 20-60 keV INTEGRAL/ISGS1 light curve of IGR J 1 6320-475 1 . 
Note that the two data sets have different lengths and the light curves have different time resolutions. The time bin of the XMM-Newton light curve is 100 s, 
while that of INTEGRAL/1SGW is about 3200 s (length of a SCW). Upper limits are given at the 3-cr level. 
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Figure 3. Lomb-Scargle periodogram of the 2—12 keV PN light curve. The 
insert shows the pulse profile of the 2-12 keV light curve obtained by fold- 
ing the PN light curve at a period of 1303 s. 
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the pulsation we restricted our period search to the 2-12 keV range 
for all EPIC cameras. The Lomb-Scargle of the PN 2-12 keV light 
curve is reported on Fig. [3] 

A very prominent peak is visible with both PN and MOS de- 
tectors. Fitting the peak with Gaussian profiles lead to best val- 
ues of 1303.8+0.9 s (7.669 ± 0.005 x 10"* Hz) and 1302.0+1.1 s 
(7.680 ± 0.006 x 10"* Hz) for MOS and PN respectively. The er- 
rors are calculated from the periodograms using the method devel- 
oped by Home & Baliunas (1986). These values are in complete 
agreement with the best value reported by Lutovinov et al. (2005), 
therefore confirming the identification of the pulsation at ~ 1300 s.A 
fainter peak, is found at a period half that of the main feature (Fig. 



Figure 4. Energy dependence of the pulse fraction, obtained after folding 
the XMM-Newton and INTEGRAL light curves. The upper limit at low en- 
ergy is given at the 3-cr level. 
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[3), that identifies it as a first harmonic of the pulse period. We then 
folded the PN light curve with a period of 1303 s (the mean of the 
of the PN and MOS values). The folded light curve is shown in Fig. 
[5] (insert). The pulse fraction (defined as P = , where l max 

and /„„■„ represent respectively the intensities at the maximum and 
the minimum of the pulse profile) is 18.11 ± 0.7% between 2 and 
12 keV for the main pulse. The pulsation is also visible in the IS- 
GRI 20-60 keV light curve. Folding the light curve at a period of 
1303 s leads to a pulse fraction between 20 and 60 keV of 17 ± 4%. 
We produced light curves in several energy ranges and folded them 
with a period of 1303 s, in order to study the energetic dependence 
of the pulse fraction. As can be seen on Fig. [4] apart from the pos- 
sible non-detection of pulsations below 2 keV, the pulse amplitude 
is compatible with being flat from 2 to 60 keV. The reported upper 
limit below 2 keV is 10.5%. 
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3.3 Spectral Analysis 

In order to avoid mixing different spectral properties together, we 
separated the observation in two different time regions, the first 
corresponding to the first flare (where we isolated ~ 3.6 ks), and 
the second corresponding to the end of the simultaneous coverage 
by XMM-Newton and INTEGRAL, i.e. the last ~ 30 ks of the 
XMM-Newton observation. The spectra were fitted between 0.6 
and 12 keV {XMM-Newton), and 20 and 80 keV (INTEGRAL). No 
relative normalization constant was included in any of the fits. 



3.3.1 Spectral analysis of the non-flaring period 

Several models were tested in the course of the analysis starting 
first with phenomenological models. A single absorbed power law 
leads to a poor reduced x~ (in the remaining of the paper xl stands 
for reduced^ 2 ). In fact a large deviation is seen at high energy, in- 
dicative of a high energy cut-off. Replacing the single (absorbed) 
power law by a power law with a high energy cut-off improves the 
fit, although the xl is still poor, and high residuals are found around 
6.4 keV. We obtain a relatively acceptable fit with a model consist- 
ing of an absorbed power law with a high energy cut-off, and a 
Gaussian emission line at ~ 6.4 keV. The xl ls 1-48 for 441 de- 
grees of freedom (dof). This value and some remaining significant 
residuals indicate that an improvement is achievable. This is partic- 
ularly true below 2 keV, where a significant excess is detected, and 
around 7 keV. To account for the latter we included an iron edge 
in the model. The fit is significantly improved with^J=1.20 (for 
439 dof). The best fit parameters are reported in Table|5|(all errors 
on the spectral parameters are given at the 90% confidence level). 
The 2-10 keV unabsorbed flux is 9.22x10"" erg cirT 2 s , and the 
20-100 keV flux is 2.33x10-'° erg cm" 2 s" 1 . 

The large excess still visible below 2 keV (Fig. [5} seems to 
be reminiscent of several of the so-called highly-absorbed sources 
detected by INTEGRAL (e.g. IGR J17252-3616, Zurita et al. 2005, 
IGR J16393-4643, Bodaghee et al. 2005), and has been seen in 
other X-ray Binaries containing pulsars (e.g. Hickox, Narayan & 
Kallman 2004). In order to try to understand its origin, we fitted 
the spectra with different types of absorption. The first one (Na ex i) 
corresponds to the absorption on the line of sight by the Galaxy and 
was modeled with the phabs model in XSPEC, while the second 
(Nh) is modeled by photo-electric absorption with variable abun- 
dance cross-sections. All the abundances where frozen to the solar 
values except that of iron given the presence of an emission line and 
an absorption edge in the simple model presented earlier. The soft 
excess is modeled by a blackbody emission. The spectral model is 
phabs* (bbody+vphabs*highecut* (power law+gauss) ) in the 
XSPEC terminology. The value of the interstellar absorption (Njiext) 
was fixed to the value obtained from Dickey & Lockman (1990) 
in the direction of the source, i.e. Nft ex ,=2. lxlO 22 cm" 2 while the 
abundances of elements are set to solar values. The fit is good with 
^ 2 ,=1.13 (438 dof). The intrinsic absorption is 1 1.8+°| x 10 22 cm" 2 , 
the blackbody temperature is 0.07^JJj keV, and the iron abun- 
dance is Z Ff ,=1.50^jj 2 j Z Q . The other parameters are compatible 
with those returned from the previous fit. The (extrapolated) 0.01- 
10 keV unabsorbed flux of the blackbody (soft excess) is 3.7x10"'° 
erg cm" 2 s" 1 , while the 0.01-100 keV unabsorbed flux of the power 
law component is ~ 4.7 x 10" 10 erg cm" 2 s" 1 . 

Since a cut-off power law is usually interpreted as a signature 
for thermal Comptonisation, we replaced the phenomenological 
model by a more physical model of Comptonisation (comptt, 



Titarchuk 1994). The choice of this model is also dictated by the 
fact that it gave a good represen tation of the XMM-Ne wton spec- 
trum of a previous observation iRodriguez et al. 20031) . although 
the lack of high energy data had prevented us from obtaining a 
good constraint on the temperature of the electrons. A simple 
absorbed comptt model (plus a Gaussian) represents the data 
rather well (^J=1.54 for 442 dof), although the soft excess and an 
iron edge are here again clearly visible. We therefore fitted the data 
with a model consisting of an (externally-) absorbed blackbody 
plus intrinsically-absorbed Comptonisation (all abundances are 
set to solar values except the iron density of the local (vphabs) 
absorption that was left as a free parameter) and a Gaussian line 
(phabs (bbody+vphabs(comptt+gauss)) in XSPEC). Similarly 
to the previous case, JVh h , is frozen to the Galactic value along the 
line of sight. We obtain a good representation of the spectrum with 
Xy=1.15 for 438 dof. The best fit parameters are reported in Table 
[3] while the "v-F v " spectrum is presented in Fig.|6| The blackbody 
(unabsorbed) bolometric flux is 2.2 xlO" 10 erg cm" 2 s" 1 , while the 
0.01-100 keV unabsorbed (i.e. external plus intrinsic absorption 
corrected) comptt flux is 4.4 xlO" 10 erg cm" 2 s . 

We also tried an alternative model involving partial covering 
by an ionised absorber (pcfabs in XSPEC). When leaving the disc 
temperature free to vary, non-realistic values are obtained for its 
normalization. We then froze kTtb to the best value obtained with 
the other models, i.e. 0.075 keV. The fit is good with xl = 1-17 
for 444 dof. The spectral parameters are compatible with those ob- 
tained with the previous model except that the intrinsic absorption 
is slightly higher here (Table|4j. Note that the fit indicates that the 
central source is almost completely covered by the absorber. 

3.3.2 The initial flare 

We started to fit the spectra with the best (first) phenomenological 
model obtained in the previous case, namely an absorbed power 
law with a high energy cut-off, a Gaussian line and an iron 
edge. The xl 15 1-06 for 439 dof. The best fit parameters are 
reported in Table [2] while the spectra are shown in Fig. [5] with 
those of the non-flare period. The 2-10 keV unabsorbed flux is 
2.39x10"'° erg cm" 2 s" 1 , and the 20-100 keV flux is 4.20x10"'° 
erg cm" 2 s"'. It is interesting to note that no soft excess is obvious 
here, while the value of _/V H has significantly decreased. It is 
possible that the soft excess is still present, although at a level 
compatible with the emission from the source (Fig.[5J. This could 
then affect the results of the fits, in particular the value of the 
absorption. In order to compare with the non-flare period, we 
fitted the data with the same types of models. We first started with 
the phabs* (bbody+vphabs*highecut* (powerlaw+gauss) ) 
model, fixing Nn ext to the Galactic value along the line of sight. 
Given the presence of an iron edge in the simple previous fit, we 
also left the iron abundance free to vary. As one could expect, the 
parameters of the blackbody are poorly constrained, and the values 
of the parameters are close to those found without the inclusion of 
the blackbody. In a second run, we fixed the blackbody temperature 
to the value found during the non-flare period (0.07 keV, with the 
phenomenological model). A good fit is achieved with ^=1.06 
(439 dof). The value of the intrinsic absorption is still lower than 
in the non-flare case with A^ H =7.8^JJg x 10 22 cm" 2 . The cut-off 
energy is now 16.2+j , keV, while the folding energy is 11. 8+} 2 
keV. The iron abundance is quite higher (although compatible 
within the errors) than in the non-flare case with Zf>=1.97^S 
The other parameters are compatible with those obtained with the 
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Figure 5. XMM-Newton/PH and INTEGRALflSGBl count spectra left: non-flare spectra right: flare spectra. In both cases the continuous line represents the 
best fit model (edge*wabs*highecut* (powerlaw+gauss) see the text for the details of the fitting). Note that the same vertical scale is used for both. 
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Table 2. Best fit parameters obtained from the fit to the 0.6-12 keV and 20-80 keV XMM-Newton and INTEGRAL flare and non-flare spectra with the simple 
phenomenological model consisting of an absorbed power law with a high energy cut-off a Gaussian and an iron edge. Errors and limits are given at the 90% 
confidence level. 
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simple model. We also allowed the temperature of the blackbody 
to vary while fixing the value of the intrinsic absorption to that 
found during the non-flare period (11.8xl0 22 cm -2 ). In this case, 
however, non credible parameters are obtained for the blackbody. 

Replacing the phenomenological model by a comptt also 
leads to a good representation of the spectra with ^J = 1.07 for 
438 dof. In a first fit the blackbody parameters are left free to 
vary. We obtain an upper limit on the temperature of 0.11 keV 
(90% confidence), but the normalization is poorly constrained. We 
therefore froze the temperature to the best value obtained during 
the non-flare period (0.075 keV) and refitted the spectra. The xl 
is 1 .07 for 439 dof. The best fit parameters are almost identical to 
those obtained when kTj,j, is free to vary, and are reported in Table 
Contrary to the phenomenological case, the N H value obtained, 
is marginally compatible with the value obtained from the fit to the 
non flare period. 

As in the non-flare case, we replaced the previous model by a 
model involving partial covering. Here again, to fit the soft excess 
well, a blackbody component is required. Since no constraints are 
obtained on the blackbody parameters if the blackbody temperature 
is left free to vary, we froze the latter to 0.075 keV. The fit is good 
with xl= 1.11 for 442 dof, and the best fit parameters are reported 
in Table0 

3.4 Phase Resolved Spectroscopy 

To study the dependence of the spectrum along the phase, we ex- 
tracted PN spectra from the non-flare period in the low part of the 
phase diagram (hereafter low phase) and from the high part of the 
phase diagram (hereafter high phase) as illustrated in Fig.Q i.e. the 
high phase corresponds to the time when the count rate is higher 
than 3.4 cts/s, while the low phase corresponds to the time when 
the count rate is lower than 3 cts/s. The origin of the phase is taken 
here as the start of the non-flare period, which we consider as the 



start of the good time intervals for the phase resolved spectroscopy. 
The shift of the phase diagram of Fig. Q as compared to that of 
Fig. is simply due to the fact that the origin of the phase in the 
former is equal to the beginning of the entire observation. Also in 
Fig-Q tne vertical axis is "real" mean counts per second rather than 
normalized ones as it is the case in Fig. [3] 

We focus on the 2-12 keV energy range, since the 0.6-2 keV 
does not show the presence of the pulsation with a rather con- 
straining 3-cr upper limit of 10.5 %. We fitted the spectra directly 
with the phabs(vphabs(comptt+gauss)) model. The Galactic 
absorption (Nn ext ) was again fixed to 2.1 x 10 22 ctrr 2 . Since we 
only consider the XMM-Newton spectra here, no constraint can be 
obtained for the temperature of the Comptonising electrons. We 
therefore fixed kT e to the value obtained during the fit to the en- 
tire non-flare spectra, i.e. 7.8 keV. In both cases the fits are good 
with^ 2 =1.18 and^ 2 = 1.01 (422 dof) for the low and high phases, 
respectively. The best fit parameters are reported in Table|5| 

It is obvious from Table|5|that the spectral parameters of both 
the high and low phase are completely compatible. The only change 
is clearly related to a change of the source flux. Although the errors 
on the line equivalent width may indicate this parameter does not 
evolve with the pulse, it is interesting to note the possible decrease 
of this parameter by ~ 30% between the low and high phase. This 
may suggest that the iron is unpulsed. 



4 DISCUSSION 

We report here the analysis of simultaneous XMM-Newton and IN- 
TEGRAL observations of the enigmatic source IGR J16320-4751. 
We focus on the time of strict simultaneous coverage by 
both satellites. We detect very significant X-ray pulsations 
at a period of aroun d 1300 s confirming previous findings 
iLutovinov et al. 2 005). The pulsation is seen in the INTE- 
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Table 3. Best fit parameters obtained from the fit to the 0.6-12 keV and 20-80 keV XMM-Newton and INTEGRAL flare and non-flare spectra with the physical 
models involving thermal Comptonisation. N\\ is the absorption intrinsic to the source, kTy, represents the temperature of the soft excess, kT,„j the temperature 
of the injected photon for Comptonisation and kT c the temperature of the Comptonising cloud. Errors represent the 90% confidence level. 
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Figure 6. XMM-Newton/PN and INTEGRAL/1SGK1 unfolded spectra left: non-flare spectra right: flare spectra. In both cases the continuous line represents 
the best fit model (phabs* (bbody+vphabs* (comptt+gauss) ) see the text for the details of the fitting). Individual components are also represented. Note 
that the same vertical scale is used for both. 
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GRALflSGRl light curve above 20 keV. Apart from the non- 
detection of the pulsation below 2 keV, no particular dependence 
of the pulse amplitude with the energy is seen. When studying the 
phase dependent XMM-Newton spectra of the source (in the non- 
flare) period, we observe no particular spectral differences between 
the high-phase and low-phase spectra. In particular, the intrinsic 
absorption, temperature of seed photons for Comptonisation, and 
plasma optical depth remain relatively constant. This is compatible 
with a model of polar accretion by a pulsar. The modulation of the 
X-ray flux is due to the misalignment of the pulsar spin axis and 
the pulsar magnetic axis. When the pulsar magnetic axis points to- 
wards us, the X-ray flux we detect is enhanced. The pulse period 
would then be the spin period as already suggested by Lutovinov 
et al. (2005). The weak amplitude of the pulse (as compared to 
other absorbed source as e.g. IGR J16393-4643, Bodaghee et al. 
2005) may indicate that the spin axis and the magnetic axis are not 
highly misaligned, or simply that the angle of the pulsar equator 
with the line of sight is quite low. The constancy of the pulse am- 
plitude above 20 keV indicates that the Comptonised component is 
also pulsed. This shows that Comptonisation occurs very close to 
the site of production of the soft photons. 

Using two independent XMM-Newton observations, we could 
refine the X-ray position to the source. This allowed us to fur- 
ther suggest that the infrared source labeled #1 in Rodriguez et 
al. (2003), and in Fig. □ here was the most likely counterpart to 
IGR J16320-4751. As discussed in Rodriguez et al. (2003), the 
magnitudes of this object could be indicative of some infrared ex- 
cess possibly due to circumstellar matter, either a hot plasma or 
some dust. 

We provide for the first time a spectral analysis of 
IGR J 16320-4751 up to 80 keV, with significant detection by IN- 
TEGRAL/ISGR1. The hard X-ray spectrum of the source is indica- 
tive of a high energy cut-off, whose parameters can not be con- 



strained when focusing on the soft X-rays only (Rodriguez et al. 
2003). The main emission mechanism in the source seems to be 
Compton up-scattering of soft X-ray photons by a hotter plasma. 
The spectral parameters we obtained either with phenomenological 
fits (Table |5J or more physical fits to the data are similar to those 
obtained for several High Mass X-ray Binaries (HMXB) (e.g. 4U 
1700-37, Boroson et al. 2003). If we assume the source is associ- 
ated with the Norma arm located between 5 and 10 kpc from the 
Sun, with a tangent at 8 kpc, we can estimate its bolometric lu- 
minosity. During the non-flare period, the extrapolated bolometric 
flux of the source is 4.12 x 10 -10 erg cirr 2 s _1 . This leads to a lumi- 
nosity of 3.15x 10 36 erg/s at 8 kpc, and 1.23x 10 36 (resp. 4.92x 10 36 ) 
for a distance of 5 (resp. 10) kpc. We note that this estimate of the 
luminosity is compatible with the typical ionising luminosity of ac- 
cretion driven X-ray (1.2x 10 36 erg/s) pulsars suggested by Bildsten 
etal. (1997). 

Our spectral analysis allowed us to reveal the presence of a 
soft excess that may be indicative of an additional medium. This 
soft excess is well fitted with a blackbody which temperature is 
quite low and found around 0.07 keV. Su ch a soft excess see ms to 
be commonly observed in X-ray binaries iHickox et al. 20041). and 
some other highly absorbed INTEGRAL s ources IZurita et al. 20051 
iBodaghee et al. 2005tlWalter et al. 20051) . This soft excess can have 
different origins from one system to another iHickox et al. 20041) . 
It could be for example the signature of thermal reprocessing of the 
harder X-rays produced by the accretion of material onto the pulsar 
either by the inner boundary of an optically thick accretion disc, or 
by an optically thin diffuse cloud. Whatever is the exact origin of 
the medium responsible for the soft excess, the fact that the tem- 
perature of the soft excess is different than that of the seed photon 
for Comptonisation rules out a model where this medium would 
be contained in the Comptonising cloud. It seems difficult to con- 
sider that this medium could be the Comptonising medium giving 
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Table 4. Best fit parameters from the thermal Comptonisation model, modified by a partial-covering absorber applied to the 0.6-12 keV and 20-80 keV 
XMM-Newton and INTEGRAL flare and non-flare spectra. Errors and limits denote the 90% confidence level. In both cases the blackbody temperature of the 
soft excess is frozen to 0.075 keV. 
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Table 5. Best fit parameters obtained from the fit to the 2-12 keV high and low phase PN spectra, with the phabs(vphabs(comptt+gauss)) model. Errors 
are given at the 90% level. 
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Figure 7. Phase diagram of the non-flare period. Only one cycle is repre- 
sented. The horizontal lines represent the limit above (respectively below) 
which the high phase (resp. low phase) PN spectrum has been extracted. 
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birth to the hard X-ray emission. Given that the Comptonised com- 
ponent shows the pulsation with an amplitude similar to that of the 
soft X-rays, it is very likely that Comptonisation occurs in the close 
vicinity of the pulsar, the seed photons being then probably emit- 
ted by the hot surface of the compact object, and the up-scattered 
pulsed Comptonised component emitted by the channeled accre- 
tion flow. 

Given the high quality of our data and the long expo- 
sure time, we detect for the first time a narrow iron line in 
IGR J 16320-4751. This kind of feature is again reminiscent 
in many HMXB, and is also detected in most of the highly- 
absorbed sources unv eiled by INTEGRAL, the extreme case being 
IGR J16318-4848 jMatt & Guainazzi 2003L IWalter et al. 2003L 
IWalter et al. 2004 iKuulkers 20051) . The presence of a narrow flu- 
orescence line, as in e.g. Vela X-l, GX 301-2, 4U 1700-37, or 
even Cen X-3 (e.g Boroson et al. 2003; Wojdowski et al. 2003), 
is usually interpreted as fluorescence of iron in a wind or circum- 
stellar matter. This indicates that accretion occurs (at least partly) 
through a wind. This points towards a high mass companion in 



IGR J 16320-4751 rather than a low mass one. The value of the 
centroid obtained during the fit to the non-flare spectra corresponds 
to Fe XIII (House 1969). We can estimate the distance R between 
the irradiating source and the inner radius of the fluorescent shell 
using £ = L/(p x R 2 ), with L the luminosity of the source, p the 
density of the gas, and £ the ionisation parameter. From Fig. 5 of 
Kallman et al. (2004) we can estimate log(£) < 2 therefore £ < 100. 
Knowing that p xR =Nh, we obtain 

R > = 10" cm ~ 0.07 AU 

N H x£ 

using the value of the luminosity found assuming a distance of 5 
kpc i.e. 1.23 x 10 36 erg/s, and R > 0.18 AU (resp 0.28 AU) for 
a distance of 8 (resp. 10) kpc. The inner edge of the fluorescent 
shell seems therefore quite far from the inner accretion flow 
which is again compatible with the presence of an additional 
medium responsible for the soft excess. It is interesting to note 
that when leaving the iron abundance free to vary, it tends to 
values greater than the solar abundance (using the values of the 
abundance of Anders & Grevesse 1989). This could be indicative 
of the iron origin and the cloud itself since both seem linked. Iron 
could have been produced by the evolution of the proto pulsar in 
IGR J16320-4751. 

Given the fact that in all the energy ranges the pulse fraction is 
roughly constant, it is very probable that the non-detection of the 
pulsation below 2 keV indicates that the X-ray flux below 2 keV, 
and consequently th e soft excess, do not pulsate. As discussed in 
iHickox et al. 2004 . this clearly rules out the accretion column as 
the origin of the soft excess. This property would suggest that the 
soft excess is either the emission by a collisionally energised cloud 
or simply the reprocessing of the hard X-rays by a diffuse cloud. 
In the latter case it seems quite natural to think that the absorbing 
material is the cloud itself. Because of energy conservation, 
F bb < AF compn where F bb is the unabsorbed bolometric flux of 
the blackbody and AF compn the difference of the unabsorbed to 
the absorbed (~bolometric) flux of the Comptonised component. 
In the non-flare case (when the blackbody components are well 
constrained), our fits lead to AF comp „ ~ 0.4 x 10'° erg cnT 2 s -1 . 
The unabsorbed bolometric flux of the blackbody component is 2.2 
xl0~ 10 erg cnT 2 s . Hence, the origin of the black body emission 
cannot be due to reprocessing of the hard X-rays. The soft excess 
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Figure 8. Error contour for the optical depth r vs. the electron temperature 
kT e from the fit with the comptt model and the absorption with variable 
abundances. The crosses mark the location of the best values, while the 
contour are the 68, 90 and 99% confidence region. 



variations of the local wind density are responsible for both the 
variations of the column density and of the hard X-ray flux. Clearly 
more observations on this and other highly absorbed sources will 
shed light on the origin of the soft excess in accreting X-ray pulsars. 
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would rather be the signature of a collisionally energised cloud. 

We also studied the spectral properties of the source during 
a flare period and compared it to the non-flare period. Although a 
strict comparison is rendered delicate by the presence of the soft 
excess which is very poorly constrained during the flare, it seems 
that the change during the flare is accompanied by a slight decrease 
of the absorption column density, while the injected temperature 
of the photons for Comptonisation and that of the Comptonising 
electrons decrease significantly. Given the degeneracy of those two 
parameters, the contour plots of t vs kT c are reported in Fig. [8] 
for both periods. This plot shows that the evolution of the Comp- 
tonising plasma optical depth and temperature is genuine. The 
ratio of the y(oc kT e T 2 , e.g. Titarchuk 1994) parameters between 
the flare and non flare spectra is y no fiare/yfiare = 0.30, indicating a 
much more efficient Compton up-scattering in the case of the flare 
period. The short time scale on which the flare occurs is more 
compatible with free-falling phenomena, rather than phenomena 
occuring in an accretion disc where times are of the order of the 
viscous time scale. This argument again points to radial accretion 
in IGR J 16320-4751, therefore further suggesting that the system 
is a HMXB. 

We note that in order to test other possibilities and explain 
the soft excess, a model involving partial covering was tested in 
the course of the analysis. However although it can lead to good 
results, it has to be noted that this model requires the addition of 
the blackbody to account for the soft excess. The result of the fits 
indicate that the covering of the central source of X-ray is quite 
high and close in both cases to its maximum value. The high values 
of the covering fraction would tend to confirm a picture in which 
the X-ray source is embedded in a dense material responsible for 
the soft X-ray. The fact that this covering fraction undergoes few 
variations between the flare and the non-flare periods seem to 
suggest that the flare is not caused by a huge geometrical change 
of this medium. This does not exclude that the flare is powered 
by inhomogeneities in the density of the stellar wind close to the 
neutron star. As the size of the region responsible for the absorption 
and fluorescence is of the order of the orbital radius, this is where 
such inhomogeneities in the wind are in fact expected. The evolu- 
tion of the absorbing column density in the model involving partial 
covering would tend to be in agreement with a picture, where 
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